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Two-dimensional 2-band insulators breaking time reversal symmetry can present topological 
phases indexed by a topological invariant called the Chern number. We propose an efficient pro- 
cedure to determine this topological index, which makes possible to conceive 2-band, tight-binding 
Hamiltonians with arbitrary Chern numbers. The technique is illustrated by a step by step con- 
struction of a model exhibiting five topological phases indexed by Chern numbers {0,±1 ± 2}. On 
a finite cylindrical geometry, this insulator possesses up to two edge states which are characterized 
analytically. The model can be combined with its time reversal copy to form a quantum spin Hall 
insulator. It is shown that edge states in the latter can be destroyed by a time reversal invariant 
one-particle perturbation if the Chern number equals ±2. 



I. INTRODUCTION 

The Chern insulators were first theoretically intro- 
duced by Haldane in 1988pJ when it was proved that 
it is possible to conceive a two-dimensional (2D) band 
model supporting integer Hall effect without a net mag- 
netic field. That was an insulator realizing the quantum 
anomalous Hall effect (QAH), with only local magnetic 
fields, breaking time reversal symmetry (TRS), and with 
bands characterized by a topological index, the Chern 
number. The presence of phases with a non-zero Chern 
number is related to the appearance of a quantified con- 
ductance at the edges of the insulator. Although a close 
relative of the regular integer quantum Hall effect, which 
is also indexed by a similar topological numberp QAH 
effect proved to be quite elusive. Only recent theoretical 
and experimental work suggests that QAH effect could 
be realized in single layer and bilayer grapheneP 

The idea of a Chern insulator laid dormant, until it 
was realized that by combining a Chern insulator with its 
time reversal copy produces a quantum spin Hall (QSH) 
insulatorPl It was proposed that spin-orbit interaction 
in graphene gaps the system such that effectively two- 
copies of Haldane's Chern insulator are realized, one for 
a given spin orientation and its time reversed mirror copy 
for the opposite orientation.^ As a result, when this sys- 
tem is put on a Hall bar, two distinct robust counter- 
propagating spin-polarized edge states show up. These 
time reversal invariant QSH insulators are characterized 
by a new Z 2 topological number P While the Chern num- 
ber is an arbitrary integer, this new topological num- 
ber indexes only two topological phases, trivial and non- 
trivial. Although the spin-orbit interaction in graphene 
was too weak to protect the edge states, soon after, the 
QSH effect based on the same mechanism has been pre- 
dicted in HgTe/CdTe quantum wella^and later confirmed 
experimentally.^ Time reversal invariant (TRI) topolog- 
ical insulators were also predicted in 3D^ and soon re- 
ported experimentally.^ A complete classification of topo- 
logical phases of one-particle fermionic insulating Hamil- 
tonians in any dimension has been achieved™ For com- 
prehensive reviews of this developing field, we point the 



interested reader to recent reviews.^ 



In the present paper, we focus on the building blocks of 
the 2D 4-band QSH insulator, namely the 2-band Chern 
insulator. We first want to point out that not all QSH 
insulators can be decomposed into two copies of Chern 
insulators. Such decomposition is only possible when the 
conservation of the third spin component s 3 is ensured. 
In such a case, the Z 2 invariant is a spin Chern num- 
ber defined by the difference of the topological indices of 
the two Chern insulators formed by each spin species.^ 
Based on geometrical considerations, we first provide a 
simple analytical method and formula to compute the 
topological index of a 2-band Chern insulator. Then we 
apply this approach to engineer general Hamiltonians for 
2-band Chern insulators. We show in particular how to 
create models with an arbitrary large Chern number, Ch, 
with only two bands. In comparison to the original Hal- 
dane modelpJ a Chern number is considered as large when 
|Ch| > 2P 



The paper is organized as follows. In Sec. |II A[ we 
describe a general procedure to analytically calculate 
Chern numbers for 2-band insulators based on the no- 
tion of degree of a map and test it on some well-known 
Chern insulators. In Sec. |III| we demonstrate how high 
Chern insulators can be engineered. In particular a new 
topological insulator with five available phases is con- 
structed. Sec.|lV]is devoted to edge states, appearing in 
the presence of boundaries. Transitions between differ- 
ent topological phases are reflected in the structure of the 
edge states on a finite cylindrical geometry. This issue 
is explored both numerically and analytically. Finally, 
in Sec. |Vj we discuss the realization of Z2-insulator by 
"doubling" the proposed model. We display the type of 
one-particle TRI perturbation that gaps the edges of the 
system composed of models in Chern number ±2 phases. 
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II. EFFICIENT DETERMINATION OF THE 
TOPOLOGICAL INDEX FOR 2-BAND CHERN 
INSULATORS 

A. Chern number as a finite sum 

A gapped 2D system with TRS breaking can be char- 
acterized by the Chern number of the occupied bands 
given bjW 

Ch=-^/ Tr(dP A PdP), (1) 

where P is the projector on the occupied bands. For 
a translationally invariant system, the Chern number is 
a topological index characterizing the U(l) vector bun- 
dle over the first Brillouin zone (BZ). A non-zero Chern 
number can be understood as an obstruction for a global 
gauge choice for the wave function on the BZ. Physically 
this phenomenon is directly related to the quantization 
of the Hall conductivity er# = Ch x e 2 /h, w here Ch cor- 
responds to the number of edge states J2H5I 

Any 2-band Hamiltonian can be expanded in term of 
pseudo-spin Pauli matrices <7j such that 

3 

J7(k) = £>(k) iff< , (2) 

z=0 

where <Jq is the identity matrix. The term h ao just shifts 
identically the energy bands and does not modify the 
topology of the Hamiltonian. We should therefore neglect 
it and write the Hamiltonian as H(k) = h(k) • er. This 
form reminds that of a Zeeman Hamiltonian for a spin 
1/2 particle in a magnetic field, except that the "field" 
h(k) is defined in momentum space. 

In the case the 2-band model the projector on the oc- 
cupied band reads P = (cr — h • <r)/2, where h = h/|h|. 
Then from Eq. ([I]), it follows that Chern number for the 
occupied band reads 

Ch=±- J^d 2 kh-(d k hxd k h). (3) 

The Chern number is a winding number that counts 
how many times does the surface traced by h wrap 
around the origin (0, 0, 0)P^ The only practical difficulty 
in determining the Chern number lies in performing the 
integration in Eq. pj). We propose instead to calculate 
Ch using a discrete summation by determining directly 
the Brouwer degree of the map h! 17 * 18 ! 

Notice that because the vector h describes an insula- 
tor with a gap at zero energy, h(k) is non-zero for all k 
in BZ. Then h(k x ,k y ) — {hi,h-i,h^) will trace a closed 
parametric surface, such that h : T 2 — > M 3 \{(0, 0, 0)} is 
a mapping from the Brillouin zone torus T 2 to a two di- 
mensional surface Ai traced by the vector h(T 2 ). Then 
h can be thought of as a composition h = tt o h, where 
7r : M 3 \{(0, 0, 0)} -> S 2 is the central projection to the 
unit sphere. 



The condition to calculate the degree of h are met: 
T 2 and S 2 are oricntable manifolds without boundary 
and have the same dimension, T 2 is compact and S 2 is 
connected. In the general case, for a point k in T 2 one 
defines the derivative map between tangent vector spaces 
dh(k) : T k T 2 -> T y S 2 . Let sgndh(k) stand for the sign 
of the corresponding Jacobian at k. Then Chern number 
is equal to the Brouwer degree 

Ch- Yl sgndii(k), (4) 

keh- 1 (z) 

where z is a regular point on the unit sphere. 

Chern number can also be computed in terms of h in 
the following way. Consider the set Y = M. H tt^ 1 (z). 
Then one has 

Ch=^ Y, sgn[(d k hxd k h)-n], (5) 

ySY keh-!(y) 

where the n is the unit vector towards z. This expression 
is just the generalization of the calculation of the winding 
n umbe r for a closed curve in 2D wrapping around a point 

The formula in Eq. ([5| can be further simplified by an 
appropriate choice of the point z. The central projec- 
tion 7r maps any intersection point between M. and a ray 
originating at (0, 0, 0) to z. If this ray does not cross the 
surface M. traced by h, it follows immediately that the 
surface does not wrap around the origin and the Chern 
number is zero. For a point z on S 2 one can immediately 
obtain a set of points on M. that project to z through 
tt. Since the expression Q does not depend on z, the 
choice of the latter can be guided by convenience. For 
instance, one can consider z lying at a coordinate axis. 
Let us choose for example the <73-axis which intersects 
A4 in a set of points. This is equivalent to say that the 
components on the other axes are zero. Said differently, 
the intersection of M. with the 03-axis are images of the 
band touchings originating from the simplified Hamilto- 
nian o\h\ + (72^2- Consequently, instead of integrating 
over the entire BZ, one only needs to consider the band 
touchings of the simplified Hamiltonian. Note that if /13 
is also zero at these points, then the system is not in a 
gapped phase. 

Now it only remains to account for the orientation of 
the surface at the intersection points with er 3 -axis. This 
can be done by studying the projection of the surface nor- 
mal vector on the (73-axis (<9fc x h x d ky h) 3 . When /i 3 > 0, 
assume that the orientation is (+1) when the sign of the 
projection is positive and (—1) when the sign is negative; 
the converse is true when /13 < 0. Then finding the Chern 
number amounts to a computation of a finite sum. Since 
we consider the entire cr 3 -axis instead of a ray, the sum 
yields Chern number multiplied by two. 

Note that although the choice of (73-axis is arbitrary, 
one must assure that the Jacobians in Eq. Q are all 
non-zero. This is equivalent to say that the band touch- 
ings are Dirac points for the simplified (two Pauli matri- 
ces) Hamiltonian. In contrast, for a gapless point with 
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a quadratic dispersion, (c\ E h x dk y h) is zero as 9fe x h is 
collinear to dk y hP3 If this is the case a different point z 
should be used. 

We can generalize and summarize the above argument 
in a formula for Chern numbers describing 2-band sys- 
tems 

Ch = - sgn(<9 fci h x d ky h) .sgn(hi) . (6) 

where i is the arbitrary axis chosen in (pseudo-)spin 
space. Therefore the integral over momenta k in the 
BZ becomes a sum over k in the set of Dirac points 
T>i for Hamiltonians H[hi = 0] (where H is the orig- 
inal Hamiltonian h • er). Note that division by two is 
required because we considered the entire axis instead of 
a ray originating in (0, 0, 0). 

Eq. ^ allows for an efficient determination of the 
Chern index and will guide us in constructing Hamil- 
tonians with arbitrary Chern numbers. Eq. ^ also tells 
us that one can start with a gapless system containing 
only two a matrices with bands that have Dirac points 
(a graphene like Hamiltonian for example). Afterwards 
one needs to compute the sign of the Jacobian Jj defined 
by 

Ji = (dkh x d ky h) i: (7) 

at the Dirac points. It contains the information about 
the Berry phase gained by a particle moving around the 
Dirac points. Note again that if the gapless points do 
not have the characteristic linear dispersion, Ji is zero.^l 
The sign of Ji gives the Dirac points "chirality". Since 
the surface M. is closed, the number of Dirac points for 
the Hamiltonians with two a matrices is always even. 

To engineer a Chern insulator it is necessary to gap the 
system by adding a "mass term". As we shall see, the 
sign of the mass term at each Dirac point can be chosen 
independently, while keeping the Hamiltonian local. As 
argued above, a necessary condition to obtain topological 
phases is that such a mass term takes different signs at 
the Dirac points. For example, if the mass term does 
not change its sign, then the surface traced by h does 
not wrap around the origin and the system is trivially 
gapped (Ch = 0). Note that this also implies the general 
statement that the sum over chiralities is zero. Finally 
it must be stressed that there is nothing special about 
this third term. Any of the other terms can play the 
role of mass for simplified Hamiltonians formed by the 
remaining ones. 

Let us call a "topological charge" the product of chi- 
rality and mass sign. From Eq. ^ we see that each 
Dirac point contributes with ±1/2 topological charge to 
the Chern number. 

From the arguments above, it is clear that a prerequi- 
site to obtain a Chern number n is starting with a Hamil- 
tonian displaying at least 2n Dirac points. Then to reach 
the largest possible Chern number (here ±n), one must 
add a mass term such that topological charges do not 



annihilate. For example if the mass sign is tuned to have 
always the same (opposite) value as the chirality at its 
respective Dirac points, one obtains maximal n (minimal 
— n) Chern number. It is clear that by tuning such a 
mass term, all topological phases with a Chern number 
between —n and n are a priori possible. Note, however, 
that in general the real space realization of this term 
might require longer than the nearest neighbor hopping. 

Graphene, for example contains, two Dirac points. To 
gap it, one can add a third term to the Hamiltonian that 
changes the sign between the Dirac points. As we shall 
see below this can yield the famous Haldane modelP^ 



B. Examples 

Before constructing a high Chern topological insulator, 
we illustrate the efficiency of Eq. ([6| to compute Chern 
numbers on a couple of popular models of topological in- 
sulators: Haldane modeiP and Bernevig-Hughes-Zhang 6 
(BHZ) "spin up" Hamiltonian for the HgTe/CdTe quan- 
tum wells. 



1. Haldane model 

Let us start by considering the paradigmatic Haldane 
modelP This model is built starting form the tight- 
binding Hamiltonian for graphene. Graphene is a hexag- 
onal lattice built out of two inter-penetrating triangu- 
lar sub-lattices A and B. It contains only nearest neigh- 
bor electron hopping with the hopping integral t\. The 
Haldane model contains also a second nearest neighbor 
hopping t%, such that when the hopping is performed 
clockwise in the unit cell the electron gains a phase <f>. 
However, the overall phase on the unit cell is zero; there 
is no net magnetic flux. Let the vectors (ai,a2,a3) de- 
scribe the displacements from B atoms to nn A atoms 
and bj = ~e.yfc(aj — a^) vectors relating nnn sites. 

The Bloch Hamiltonian reads 

3 

-ff (k) = {2i 2 cos((/>) cos(k • bi)<7o + ti[cos(k • ai)<7i 

i=l 

M 

+ sin(k • ai)cr 2 ] + [y - 2t 2 sin(</>) sin(k • b|)]o" 3 }, 

(8) 

where ±M is the on-site energy. Let us study the topol- 
ogy of the surface traced by vector h. Again, if one 
chooses to study the points where the surface inter- 
sects er 3 axis, it follows immediately that those are ex- 
actly given by the Dirac points of graphene, i.e. by the 
Hamiltonian with M and t 2 — > 0. It is easier to 
reparametrize the Hamiltonian, such that k • b3 = q x 
and k • h ± = q y . This describes a mapping to a square 
BZ, with Dirac points given by q± = ±(27r/3, 27r/3). Be- 
cause SLi = §ey fe (bfc - b,-), k ■ (ai,a 2 ,a 3 ) = -\(-2q x - 
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QyiQx — Qy,<lx + 2q y ). Dirac points' chirality is given by 
the sign of the Jacobian in Eq. Q. We find 

sgn(J 3 ± )= T l, (9) 

where denotes the Jacobian calculated at the Dirac 
point q^.. 

In order to calculate the Chern number one must also 
consider the mass sign at the Dirac points, m±(q±) = 
M =F 3v / 3i 2 sin(0). Therefore, following Eq. m one re- 
covers Haldane result for the Chern number 

Cn = ^[ s g n ( m -) ~ sgn(m+)]. (10) 

The transition from a topological insulator to a nor- 
mal insulator is marked by a semi-metal state, where 
the gap closes at least at one Dirac point when M — 
±3\/3i2 sin(0). In the normal insulating phase Dirac 
points have identical associated mass sign, such that Ch 
goes to zero. 
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TABLE I. Values of the Chern number Ch according to the 
different values taken by M and B > 



It follows from Tab. U that as M varies between and 
SB, the Chern insulator H(k) exhibits two topologically 
nontrivial phases with Ch = ±1. When M is outside 
(0, 8B) region there is only a trivial insulator phase. Note 
that with respect to -ff(k), its time reversed copy will 
always have the opposite Chern number. 



2. BHZ model 

The Chern number calculation can as well be exempli- 
fied in the case of the recently discovered Z 2 insulators 
such as the 2D HgTe/CdTe quantum wells^The low en- 
ergy Bloch Hamiltonian is written in a basis of four states 
\El,mj = 1/2), \Hl,mj = 3/2), \El,mj = -1/2), 
\H1, mj = —3/2) and it has the form 



H 



(H(k) 
^ £P(-k) 



(11) 



The system is assembled out of two Chern insulators, 
H(k) and its time reversed copy. To illustrate the calcu- 
lation of the Chern numbers using a discrete sum, it is 
enough to pick one Chern insulator: 

H (k) = A s'm(k x )(Ti + A sin(fc a )(72 

+(M - 2B(2 - cos(k x ) ~ cos(k y )))a 3l (12) 

where A, B, M are parameters of the Hamiltonian. Let 
us consider again the surface traced by h and choose 03 
as a special axis. The points where the <73-axis pierces 
the surface are given by the condition that hi and hi 
terms vanish. This gives four "Dirac points" (k x ,k y ) 6 

{(0,0),(0,7r),(7r,0),(7r,7r)}. 

The chirality of each Dirac point is given by the sign 
of Jacobian J3 



sgn( J 3 ) = sgn[cos(fc :E ) cos(k y 



(13) 



evaluated at the Dirac points. The mass term /13 has the 
following expression at the Dirac points, /13(C), 0) = M, 
h 3 (0,ir) = M^O) = M - 2B and /i 3 (7r,7r) = M - AB. 

The Chern number can then be easily computed for 
different values of M and B. The results for the case 
B > are summarized in Tab. JJ 



III. HIGHER CHERN NUMBER MODELS. 
CONSTRUCTION OF A CHERN INSULATOR 
WITH PHASES {0,±1,±2} 

In the following section we present the construction of 
a 2D model with a higher Chern number. The strategy 
will be to work with an abstract Hamiltonian in fc-space. 
We delay its implementation on a particular lattice at 
the end of this section. New topological phases are made 
possible by tuning the parameters of the mass term and 
controlling the addition and subtraction of "topological 
charges" . 



A. Topological phases in fc-space Hamiltonian 

The most important common point between Haldane 
model and "half" -BHZ Chern insulators is the existence 
of only three topological phases described by Chern num- 
bers in {0, ±1} (0 characterizing the trivial phase). Using 
Eq. (6| it is simple to construct an artificial model which 
has larger Chern number. We assume that we are in 
momentum space with a general 2D Hamiltonian of the 
type 



hi{k x , ky)ai + h 2 (k x ,k y )(j2- 



(14) 



h\ and h 2 are continuous, 27r periodic functions of k x and 
k y . 

The bands meet when both functions hi and h 2 are 
zero. This defines a system of two equations with two 
variables k x ,k y , whose solutions define points in fc-space. 
As explained above, in order to obtain a higher Chern 
number (here ±2) it is necessary to have at least four 
Dirac points. 
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TABLE II. Chirality \ °f the four different Dirac points 



Let us start by considering the basic template on which 
the topological insulator is built. Among the simplest 
models, one can consider the following Hamiltonian: 

H (k x ,k y ) = 2t 1 [cos(k x )a 1 + cos(k y )a 2 }. (15) 

Here, we assume that the Pauli matrices a correspond to 
some pseudo-spin or orbital degree of freedom. Note that 
already the system breaks TRS, THq^T- 1 ^ H (-k). 
The energy dispersion reads 




E = ±2tiJcos 2 (k x ) + cos 2 (k y ), 



(16) 



such that there are four Dirac points ko = (±7r/2, ±7r/2). 

The Dirac point chirality is given by the sign of the 
Jacobian J3 



sgn(J 3 ) = sgn[sin(/c :I .)sin(fc y )]. 



(17) 



This determines immediately the chirality \ of the four 
Dirac points as summarized in Tab. Ill) 



Notice that Dirac points at k and — k have the same 
chirality. They will be referred in the following as a pair 
of Dirac points. 

To obtain a Chern insulator one needs to add a mass 
term. As we can see from Eq. the Chern index 

depends on both the chiralities of the Dirac points and 
the sign of the mass term in their vicinity. Let us add a 
mass term of the form h^{k x ^ k y )a^. Since /13 is a periodic 
function on the BZ, in the general case its zeros form a 
set of closed lines on the two-dimensional torus. The first 
condition in order to gap the initial system is that these 
lines must not pass through the Dirac points. Thus the 
lines of zeros delimit regions, TZ\ and 72.2, where the mass 
term has the same sign. For a proper choice of the gap 
term it is somewhat simpler to see the problem from a 
geometrical point of view. In order to maximize |Ch|, one 
needs that the lines of zeros separate the pairs of Dirac 
points such that a pair of points of a given chirality are 
contained in a region of positive mass, while pair of points 
of opposite chirality are contained in a negative mass 
region. In short, each pair of Dirac points are placed in 
regions where the mass term has different values. On the 
contrary, if a pair of Dirac points is "broken" (such that 
one is in region TZ\ and the other in region 7^2), then 
the topological charges will cancel out as can be directly 
inferred from Eq. Q. 

Let us first realize the topological phases with a Chern 
number ±2 where each pair of Dirac points is in dif- 
ferent regions 1Z. The periodic function cos(k x + k y ) 
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FIG. 1. (a) Domains associated to the mass term cos(k x + k y ). 
Phases with Chern number ±2 are realized, (b) Domains asso- 
ciated to the mass term sin(k x ) +sin(k y ). Phases with Chern 
number ±1 are possible after the remaining Dirac points are 
identically gapped. In both figures, dashed lines denote the 
locus of points of zero mass term, while IZi and IZ2 denote 
regions with opposite values for the mass term. The solid cir- 
cles represent Dirac points of the initial gapless system; the 
corresponding chiralities are marked by the signs within. 



accomplishes such demand with lines of zeros given by 



ky — 



271+1, 



term are represented in Fig. ( la 



The regions TZi and 72-2 for this 
The mass term is neg- 
ative in region 1Z\ and positive in region 7?-2- The lowest 
Chern number Ch = —2 is obtained for this model. Note 
that the above choice of /13 is not completely arbitrary, 
since the term i 2 cos^ + k y )a^ corresponds to a short 
range hopping with the amplitude ti- Hence, the modi- 
fied Hamiltonian H\ becomes 

i?l(k) = 2^003(^)171 + 2^1 COS (fcj,) 172 +2^2 COs(fc x + fcj,)(T3. 

For this Hamiltonian, one readily computes Ch — 
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2sgn(— t-i). Therefore the sign of t% determines whether 
a maximal or minimal Chern number is reached. Note 



that the Hamiltonian in Eq. ( 18 ) lacks trivial insulating 
phases or Chern phases with Ch = ±1. 

To produce a trivial phase with Ch — 0, it suffices 
to add a large, "staggered potential" ma 3 such that the 
entire mass term for all Dirac points has the same sign. 
Note that since Pauli matrices can act on any degree of 
freedom, the term "staggered" could for instance refer to 
the fact that two sites in a cell or two orbitals on a site 
have an associated ±m constant energy. 

The Hamiltonian then changes from Hi to H2 = Hi + 
ma 3 and consequently the Chern number becomes 



Ch = sgn(— to — 2i 2 ) + sgn(ro — 2t 2 ) 



(19) 



Therefore, when the mass is large enough, |m| > 2^1, 
the system enters a trivial phase. We have therefore a 
transition from a Ch = ±2 phase to a Ch = phase. 

Along the same reasoning, we can obtain phases with 
Ch = ±1. Remember that Dirac points at k and — k have 
the same chirality. Then to obtain a Ch = ±1 it is neces- 
sary that the mass term is allowed to take different signs 
for exactly one pair of Dirac points. If one adds only even 
functions of k to the mass term, it would not be possible 
to "break" a pair of Dirac points. Consequently, an odd 
function is required. The simplest choice would be to add 
the term proportional to sm^fe^) + sin(fc B ). The mass for 
one pair of Dirac points is unchanged, while for the pair 
(ko,— ko) with ko = (tt/2,tt/2) the mass changes (see 
Fig. lb for the model containing only sin(k x ) + sm(k y ) 
in the mass term). Note that as follows from Fig. |lb| if 
the mass term contains solely sin^) + sin(fc y ), the sys- 
tem is not an insulator because two Dirac points are not 
gapped. However, the presence of even functions in the 
mass term gaps identically the two Dirac points so that 
phases with Chern number ±1 become possible. Adding 
all the terms together gives the following complete Bloch 
Hamiltonian H 

H = 2ti cos(k x )<Ji + 2ti cos(ky)<j2 

+[m + 2ti cos(k x + k y ) + 2t3(sin(fc x ) + sin(fcy))]<73. 

(20) 



There are four free parameters (to, ti, £2, £3) in our model, 
and five insulating topological phases {0, ±1, ±2}. Let us 
assume that all parameters are real. All the phases can 
be reached by varying (to, £2,^3) while keeping ti fixed. 
The Chern number for the final model reads 



Ch = sgn(-TO - 2i 2 ) + - [sgn(TO - 2< 2 + 4* 3 ) 
+sgn(TO - 2i 2 - 4i 3 )]. 



(21) 



Such a formula gives immediately the phase diagram 
associated with the system described by the Hamiltonian 
in Eq. ( 20 ) . Notice that the parameter ti does not enter 



in the description of the phases, but manifestly needs to 
be finite to have non-vanishing <j\ and 02 components in 
the model. The formula (21) can be illustrated by the 
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FIG. 2. Phase diagram of the system for t% = 1. Each region 
is denoted by the corresponding Chern number and represents 
an insulating topological phase. The boundaries of the regions 
represent topological transitions where the system becomes 
gapless. (color online) 



phase diagram on the Fig. This diagram contains 
only four phase of the model; the phase with Ch = 2 
needs t 2 < 0. 



B. Direct space realization 

For the moment the system was assumed translation- 
ally invariant in 2D and was described abstractly in pa- 
rameter space. However to investigate the edge states it 
is necessary to propose a lattice implementation for the 
model. 

The system could be realized on a triangular lattice 
with two orbitals on each site. The parameters are inter- 
preted clS kx^y — k ' 3-1.2 

with ai and a.2 the two Bravais 
vectors. The Pauli matrices are then operating in the 
orbital space. 



The Hamiltonian ( 20 ) can be rewritten as 



X — T r f ^ t / \ 

+4 j+ i(h(7 2 + it 3 a 3 )cij + c\ +lj+1 t 2 (T 3 Cij + h.c] 



where 



Ck 



(22) 



(23) 



The position is defined by r,y = iai + ja.2 with a.1/2 
as Bravais vectors making an angle 27r/3 between them 
(see Fig ([3])). We set the lattice constant a = 1. On 
each site there are two different orbitals with energy ±to. 
The Hamiltonian ( 22 ) describes through i-terms the over- 



lap between (iniequivalent orbitals. The model is repre- 
sented in Fig. (3j). Note that there is no net flux per- 
pendicular to the two dimensional plane. But TRS is 
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(i + U) 



FIG. 3. Direct space realization of the model (20 1 
Eq. (22 1. denotes lattice sites 



sec 



11/2 



arc Bravais lat 



tice vectors. • (o) represents orbitals with energy m (— m). 
The vertical axis represents the on-site energy difference be- 
tween the two inequivalent orbitals. Black lines represent t\ 
hoppings, blue lines, hoppings, and red lines, tz hoppings. 
An arrow on a link indicates that an electron hopping in the 
corresponding direction gains a 7r/2 phase. Similarly a double 
line indicates a 7r phase gain, (color online) 



broken, because for certain closed paths the electron can 
still gain a non-zero phase. 



IV. BULK-EDGE CORRESPONDENCE 

Until now, only bulk properties of the model were ex- 
plored. In the present section, we investigate how these 
properties manifest themselves through the edge state 
structure. It was proved for general 2D Chern insula- 
tors that the Chern number is reflected in the number of 
gapless edge states one encounters on a finite geometryPH 

We have therefore implemented the aforementioned 
system on a cylinder and investigate subsequently the 
edge state wave functions and edge state dispersion at 
its two ends. In the first subsection, mainly numerical 
results are presented. It illustrates the fact that the 
Chern number reflects itself in the number of gapless 
states at the interface between phases with Ch = — 1 and 
Ch = — 2. In the second subsection, analytical methods 
are employed to obtain the edge state structure. Here 
the main goal is to obtain the edge wave function and 
dispersion for any parameter values. In particular, the 
complete edge state determination for a choice of param- 
eters in high Chern —2 phase is given. 



to a2 are glued together to obtain the the final cylin- 
drical shape. Because translational invariance is main- 
tained in the direction parallel to a! , k x remains a good 
quantum number. Therefore, one can fix k x and con- 
sider the resulting one-dimensional problem. Let us write 
the one-particle solutions of the corresponding stationary 
Schrddinger equation as 

m x ,j))=Y,Uk x )c\j$), (24) 

fex 

where j denotes the layers of sites in a2 direction, and 
ijjj is a spinor due to the fact that there are two orbitals 
in the problem. 

Then for a given quasi- momentum k x , the Schrddinger 
equation reads 



Eil>j, 



where 



Ti = ma 3 
T 2 = t\Oi ■ 



■ 2ti cos(fc x )<7i + 2t 3 sin(fc x )(73, 
(t 2 e lk * -it 3 )a 3 . 



(25) 



(26) 



The cylinder has the edges at j = 1 and j 
Hard wall boundary conditions are imposed, such that 
the amplitudes i/jq, i/'l b +i must vanish. The dispersion 
E as a function of k x is obtained by numerically solv- 



ing Eq. ( 25 ) for the given boundary conditions and for 
different choices of the parameters (see Fig. Q). In 
our numerical experiments the cylinder circumference is 
L x — 180 sites and with height L y = 40 sites. All the 
energies are measured in units of t\. 

To illustrate the bulk-edge correspondence, we sample 
several regions of the phase diagram on Fig. (|]). In non- 
trivial topological regions, edge states appear around the 
ends of the cylinder. From the phase diagram we can 
predict the number of edge states in our model. For ex- 
ample, three sets of parameter values along the constant 
m = —1.4 line are taken such that the transition between 
the —1 and —2 phases is explored (see Fig. (4d)). 

While the bulk remains insulating there are states 
crossing the gap. These are the edge states and their 
total number is 2 x |Ch| since the cylinder has two edges. 
Note that the edge states at zero energy, cross the gap 
at k x = ±7r/2. At any topological transition the bulk 
closes at least in one of the special points k x = ±ir/2. 
A transition changing the Chern number by two requires 
that the gap closes at both points, while for a change of 
one, only one Dirac cone forms. 



B. Analytical solution 



A. Numerical experiments 

The finite geometry chosen for the numerical study 
is a cylinder. It is constructed out of a patch of the 
lattice having the shape of a parallelogram Bravais vec- 
tors ai and a2 as edges. Subsequently the edges parallel 



A greater insight into the model is gained by solving 
Eq. ( 25 ) analytically. It allows one to clearly identify the 



gapless states as edge states and also to determine their 
penetration length. The edge state dispersion law can 
be found either by a direct solution of Schrddinger equa- 
tion or indirectly by the method described in Ref. |22j . 



(d) 
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FIG. 4. Energy dispersions for (a) ts — 0.4 (Ch = —2), (b) t$ = 0.85 (topological phase transition at closing bulk gap), (c) 
t$ — 1.6 (Ch = —1). The other parameters are t\ = 1, ti = 1, m = —1.4. The simulation is done when the system is on 
cylinder with height of 40 sites and a circumference of 180 sites. The number of edge states is 2 x |Ch| because there are two 
edges, (d) Representation of chosen points on the phase diagram, (color online) 



We explore both approaches in the general setting of our 
model and exemplify the results for a particular choice 
of parameters corresponding to the phase Ch = —2. 

As it was elegantly proved in Ref. [55] the condition 
of existence and edge state dispersion can, under certain 
provisions, be found from a simple analysis of bulk Bloch 
Hamiltonians. The method developed there applies when 
an infinite ribbon or a cylinder is cut out of the infinite 
2D system. The direction of the cut must follow a Bravais 
lattice vector. In this case the momentum parallel to the 
cut fcy is conserved and the system splits into a set of 
ID Hamiltonians describing the motion of the electron 
between the layers of sites parallel to the cut. The final 
prerequisite to apply the method is that there are only 
nearest neighbor layer hopping terms. Eq. (25) shows 



that it is indeed the case in the present model with k\\ = 
k x . 

According to Ref. [55], the key information about the 
edge states can be revealed by studying the curves traced 
by h as a function of fcj_ with fixed k\\ . In the case 
when only the nearest neighbor intcrlaycr hopping is al- 
lowed these curves are planar (actually they are ellipses) . 
Therefore, h can be decomposed in two parts, hj_ per- 
pendicular to the ellipse plane and hn the in-plane com- 
ponent. Each component yields some important piece 
of information about the edge states. Namely, the edge 
state with a given fey exists if and only if the ellipse traced 
by h\\ encloses the projection of the origin onto the plane 
of the ellipse. The energy of the state is equal to ±|hj_|. 

In our case — k x and k± — k y . This yields 

h|| = (0, 2ti cos(fcj,), 2t 2 cos(k x ) cos(ky) 

+2{t 3 - t 2 sin(k x )) sm(ky) +m + 2t 3 sin^)). 

(27) 

For a fixed k x , the equation ( |27| ) describes an ellipse 
parametrized by k y G [0, 27r). The condition that the 
ellipse encloses the origin reads 



\m + 2t 3 sm(fe x )| < 2|t 3 -t 2 sm(k x )\. 



(28) 



This equation determines the range in k x where edge 
states exist. The energy of the state is ±2ti cos(k x ). 



Although the edge dispersion is determined, it must 
accommodate up to four edge states. The existence con- 
dition also gives for |Ch| = 2 two intervals of allowed k x 
where edge states exist. 

The edge states can also be explored directly by study- 
ing the special solutions of equation (25). The strategy 



consists firstly in finding a zero energy solution. From the 
above analytical results or from numerical experiments it 
is apparent that if there are edge states, they will cross 
zero energy only at k x = ±7r/2. After obtaining the 
solution at this particular fc^-points, one extends the so- 
lution for the entire range of k x which allows edge states. 
Note that although this method relies on the starting in- 
formation about the zero energy crossing of edge states, 
it determines the existence condition and dispersion by 
itself. 

Let us consider the solutions to Schrodinger equation 
( 25 1 of the form ipj — p J <fi. The complex parameter p cap- 



tures the behavior of the wave function in the direction 
transverse to the edge p = e lky and 4> carries the spinor 
structure of ipj. The ensuing equations, at k x — ±n/2 
read 



(0,0) T = {[(h<J 2 -i{t 3 ^t 2 )a 3 )]p 

+ [{t 1 cr 2 +i(t 3 Tt 2 )a 3 )}p- 1 



(±2t 3 



m)a 3 }(f>. 
(29) 



By multiplying Eq. (29) with a 3 on the left hand side, 
we obtain an equation which depends only on a± . There- 
fore 4> is a linear combination of the eigenstates |a;±) of 
o i . The solution can then be extended beyond the special 
points k x = ±7r/2 by continuity. 

Eigenmodes p are found from the general fc^-dependent 
equation 



(r 1 + r 2j0 - 



x±) 



0. 



(30) 



Note that there are two solutions p for a given eigenstate 
|ir±). A certain symmetry between the solutions is im- 
mediately apparent. If p\ and p 2 are solutions to the 
equation for the eigenstate \x+), then p*~ and p 2 ~ are 
solutions for the equation corresponding to \x— ) eigen- 
state. In particular, p\ and p 2 associated to \x+), are 
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determined from 



0= TO + 2i 3 sin(fe r ) + t 2 (e lK p + e 



i k j 



p- 1 ) - it 3 ( P - p- 1 ) 
(31) 



Note that the edge state dispersion follows as a byproduct 
E = ±2ti cos(k x ) depending on the eigenstate |x±). 

For a given fcn the general solution to Eq. (251 thus 
reads 



i>j = (ciPi + c 2 pi)\x+) + (c 3 pl 3 + c 4 p* 2 ~ 



(32) 



This result is further constrained by imposing the bound- 
ary conditions. In order to find localized wave func- 
tion near j = 1, one imposes the boundary conditions 
4> = (0, 0) T , such that the general solution now writes 

^ = c + (p{ 4)\x+) + c-(pi-* - pt j )\x-). (33) 

For an edge state, one cannot have \p\ — 1 since this 
would correspond to a freely propagating mode perpen- 
dicular to the edge. The localized solutions at j = 1 exist 
if and only if one of two following conditions is satisfied: 



\pi\ > 1= 



\Pi\ < 1, 



|P2|>1, C + =0, 



o. 



(34) 



(35) 



Thus if there is an edge state solution, its form and 
energy are 



E = 2ti cos(fc x ) 



(-) 



P-2 



E = —2t\ cos(k x ). 



where 



with 



Pi, 2 



-b ± Vb 2 - Aac 
2a ' 



(36) 



(37) 



(38) 



a = —it\ + t 2 e * - it 3 , 
b = m + 2i 3 sin(fc a; ), 



c = —%t\ + t 2 e 



—ik x 



+ lt 3 - 



(39) 



larger then one and the other is smaller than one. When 
|Ch| = 1, only one of the solutions holds in the Brillouin 
zone. When |Ch| = 2, there are k x where both solutions 
hold. 

Let us illustrate the above results for a special point 
ti = l,t% = l,t3 = 0, m = of the phase diagram p]). 
This point corresponds to the "center" of the Ch = —2 
phase and is characterized by the largest gap and flattest 
bands for the spectrum of the bulk states. One can ex- 
pect two edge states at either end of the cylinder. The 
eigenmodes p are determined by 



Pi,2 = ±i 
Therefore 

|Pl,2| 



—iti + t 2 e lk * 



t\ + tl + 2txt 2 sm(k x ) 



1/4 



(40) 



(41) 



t\ + t| — 2t\t 2 s\a(k x 
and the localization length £ of the edge state is given by 
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-1/MIPi.al). 



(42) 



It is apparent, either by plotting the eigenmodes' ab- 
solute value or by checking the existence condition of 
Eq. (28 1, that edge states are expected for any k x . There- 



fore this choice of parameters proves to be quite partic- 
ular, giving maximal extension for edge states. 

Evaluation of \p\ shows that \p\ > 1 for k € [0, ir] and 
\p\ < 1 for k £ [— 7T, 0]. That means the edge state local- 
ized near j = 1 (bottom) is an eigenstate \x—) in the BZ 
interval (0,tt) and the other is an eigenstate of \x+) in 
(— 7r, 0). The j — L y solution (top) immediately follows 
by symmetry. 

As expected, there are four edge states, with wave func- 
tions and energies given by: 



-Kb 



c+ (i-(-i)Vik+) 



E = 1t\ cos{k x ) 



</f 6 = c_(l-(-l)>Pk-> 
E = —2ti cos(k x ) 

i;f t = c + (l-(-iy)pi\x+) 

E = 2ti cos(k x ), 



-A't 



c_(i-(-iy)p^\x-) 



E = —2t\ cos(k x 



k x G (-7T,0) 
k x G (0,tt) 
k x G (0,tt) 
fe€(-7r,0). (43) 



The indices t and 6 indicate whether the edge states live 



In summary, there can be at most two solutions local- 
ized at j = 1. The existence conditions (34) and (35) are close to the top (j — L y ) or the bottom (j = 1) part 

of the cylinder. ±K indicates whether the edges state 
crosses the zero energy at ±7r/2 or, equivalently in this 
case, whether it is extended in the right, respectively left, 
part of the BZ. The coefficients c± are normalization co- 
efficients which are not of interest here. The edge states' 
spectra (43) are plotted in Fig. ([5| together with the nu- 
merical solution. 



imposed to determine the number of edge states and their 
extension. Note that for every solution at edge j = 1 
there is the solution at j — L y which can be readily ob- 
tained by inversion and conjugation of p's. 

The topological properties of the bulk phase manifest 
themselves as existence of edge states at j = 1. When 
Ch = no edge state solution exists, since one of \pt\ is 
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fcq o 




FIG. 5. Energy spectrum as a function of momentum k x on 
a cylindrical geometry (height 40 sites and circumference 180 
sites). Two edge states located near j = 1 are represented in 
blue, and two at j = 40 in green. The position and chirality 
of the edges is schematically represented in the inset, (color 
online) 



the conditions for the energy gap are not changed. That 
means the previously found insulating phases remain in- 
sulating phases for the new model. 

Edge states are again found when going to the cylindri- 
cal geometry explored in the previous section. Because 
there are no spin mixing terms, the energy spectrum 
is trivially obtained by "doubling" the spectra already 
found for the spinless Hamiltonian. More precisely it is 
obtained from the union of the spinless (now spin up) 
Hamiltonian spectrum and its reflection about k x = 
under TRS. Therefore the number of edge states will 
also double such that each original edge state will get 
its Kramers partner. 

Although every previously nontrivial Chern phase will 
show edge states in the spinful model, not all of them 
are robust. Indeed, as we shall see, the QSH insulator 
constructed out of the spinless model with Ch = ±2 al- 
lows for a one -part icle TRI perturbation that destroys 
the edge statesPISH 

The low-lying edge states are described by an effective 
Hamiltonian, obtained by linearizing the solutions (43) 
near k x = ±7r/2. At a given edge, for the phase with 
Ch = ±2, this yields: 



Notice that the pair of states living at either end of the 
cylinder have the same direction of propagation. Taking 
the derivative of their dispersion shows that the velocity 
is identical. Still they are distinguished by a "valley" 
quantum number ±K, determined by where they cross 
zero energy. 



V. Z 2 INSULATOR 

In this section, a quantum spin Hall (Z2) insulator is 
created out of the Chern insulator in Eq. ( 20 ) by adding a 



spin flavor to the model. In particular, we show that the 
edge states constructed from the ±2 Chern numbers are 
not robust and actually correspond to a trivial insulator, 
in agreement with the Z2 classification. This is confirmed 
by the existence of one-particle, TRI, local perturbations 
that open the spectral gap. 

We suppose that the third component of the spin S3 
commutes with the Hamiltonian. Thus the system con- 
sists of two independent components. We suppose that 
the spin up component is described by (20), while the 
spin down one represents its time reversed copy. This 
yields the following 4-band Hamiltonian: 

H(k) = 2ti cos(k x )<Ji (g) so + 2ti cos(fc y )(72 ® S3 
+[m + 2t 2 cos(k x + k y ))a 3 <g> s 
+2i 3 (sin(fc x ) + sm(k v ))a 3 <g> s 3 , (44) 

where s represents electronic spin and ao, so are identity 
matrices. 

Because the spin Hamiltonian is created by two copies 
of the spinless Hamiltonian, with no spin mixing terms, 



H eS {q x ) = ^ Km vq x ^ K R/t + ^_ KRJ[ vq x ^^KB^ 

- ^-KLi^-KLi - & KLi vq x V KLl , (45) 

where Fermi velocity reads v — 1t\. The indices of the 
fermionic creation and annihilation operators ^ and ^> 
describe the valley (±K), the direction of motion (L or 
R) and that of spin (t or |) . Note also that the first two 
terms in H e g describe the dynamics of spin up electrons, 
and therefore correspond to the original 2-band Hamil- 
tonian, while the spin down terms stem from of time 
reversal operator T; T^ Km T~ l = ^ —kli- The lock- 
ing between the direction of the spin and that of motion 
means that iJ e ff describes a helical liquid. 23 

The edge states above are not robust because one can 
create the following one-particle, TRI, local perturbation 



that will gap the edge helical liquid in Eq. ( 45 1 (local 
meaning there is no inter-edge scattering) 



h.c. 



(46) 



It is possible to build many tight-binding perturbations 
yielding the above form at low energy. It is noteworthy 
to observe that they all break the spin S3 symmetry. An 
example of TRS perturbation in the tight-binding formu- 
lation is £4 sm{k x )a 3 ® s\. 

For the phases with Ch = ±1 no one-particle, local, 
TRI perturbation can result in backscattering of the edge 
states. The above example agrees with the statement 
that models with an even number of Kramers pairs of 
edge states are Z2-trivial. 12J 
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VI. CONCLUSION 



We have developed a way to calculate the Chern num- 
ber by a discrete summation for 2-band, two-dimensional, 
TRS breaking, topological insulators in a tight-binding 
formulation. The method allows for an efficient determi- 
nation of topological phases in such systems and helps to 
conceive models with large Chern numbers without in- 
creasing the number of energy bands. We start with sim- 
ple gapless Hamiltonian containing only two Pauli matri- 
ces and 2n Dirac points. Then the system is gapped by 
adding an explicitly designed mass term containing the 
third Pauli matrix, yielding any of 2n + 1 possible topo- 
logical phases (the Chern number can thus vary between 
{-n,-n+l,...,n}). 

This technique is illustrated by constructing an "arti- 
ficial" Chern insulator which has five available topolog- 
ical insulating phases described by the Chern numbers 
{0,±1,±2}. After determining completely the phase di- 
agram and the bulk properties of the model, we have 
studied the edge structure of the model on a finite cylin- 
drical geometry. For the particular choice of edge geome- 
try, the conditions for existence of gapless edge states and 
analytical expressions for their dispersion were obtained. 
The edge wave function structure is also determined for 
the edge states in any topologically nontrivial phase. For 



the case Ch = —2, we consider in detail the edge state 
solution and show it to be in agreement with numerical 
simulations. Finally we discuss a 4-band realization of 
a QSH Z2-insulator from the 2-band Chern insulator at 
hand. We determine the phase diagram of the 4-band 
model and show that the Chern phases ±2 indeed yield 
the edge states unstable with respect to one-particle TRI 
perturbations. 

The technique described in the present paper could 
be used to design 2-band tight-binding models featur- 
ing the QAH effect with large Chern number, realiz- 
able with ultracold fermionic atoms in an optical lattice. 
Recently, several 2D optical lattices — honeycomb, brick- 
wall, kagome, checkerboard, etc. — with more than one 
orbital per unit cell were realized experimentallyPStHl j n 
addition, methods to create artificial gauge potentials- 
and therefore complex hopping amplitudes — are now 
reaching maturity.^ Thus it appears feasible to tailor an 
optical lattice that could feature non-trivial Chern insula- 
tor. The result could be verified experimentally by mea- 
suring the Berr y cur vature through Bloch oscillations, as 
recently shownpSEH or time-of-flight measurements.^ 
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